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Abstract. We present VLT/ISAAC near-infrared spectroscopy of two gravitationally-lensed z ^ 1.9 galaxies, A2 and S2, 
located behind the cluster AC 114. Thanks to large magnification factors, we have been successful in detecting rest-frame 
optical emission lines (from [O II] A3727 to Ha-i-[N II] A6584) in star-forming galaxies 1 to 2 magnitudes fainter than in previous 
studies of Lyman break galaxies (LBGs) at z ^ 3. From the Ha luminosity, we estimate star formation rates (SFRs) of 30 and 
15 M0 yr~^ for S2 and A2 respectively These values are 7 to 15 times higher than those inferred from the UV continuum flux 
^ ■ at 1500A without dust extinction correction. In setting SFRhc^ ~ SFRuv, one derives extinction coefficients E(B-V) ^ 0.3 for 

' S2 and E(B-V) ~ 0.4 for A2. The behavior of S2 and A2 in terms of 0/H and N/0 abundance ratios are very different, and 

, they are also different from typical LBGs at z ^ 3. S2 is a low-metallicity object (Z ^ 0.03 Zq) with a low N/0 ratio, similar 

■ to those derived in the most metal-poor nearby Hll galaxies. In contrast, A2 is a high-metallicity galaxy (Z ~ 1.3 Zq) with a 

I high N/0 abundance ratio, similar to those derived in the most metal-rich starburst nucleus galaxies. The line-of- sight velocity 



o 



dispersions, derived from emission line widths, are 55 and 105 km s , yielding a virial mass of 0.5 and 2.4 xl0^° M0,for S2 
and A2 respectively. Thanks to the gravitational amplification, the line profiles of S2 are spatially resolved, leading to a velocity 
gradient of ±240 km s~^, which yields a dynamical mass of 1.3 x 10^° Mq within the inner 1 kpc radius. Combining these 
new data with the sample of LBGs at z ^ 3, including the lensed galaxy MS 1512-cB58, which is the only LBG for which 
Oh' physical properties have been determined with similar accuracy, we conclude that these three galaxies exhibit different physical 

Q , properties in terms of abundance ratios, SFRs, M/Lb and reddening. High-redshift galaxies of different luminosities could 

^ ■ thus have quite different star formation histories. 

to 



Key words, galaxies: evolution - galaxies: starburst - galax- and limited to the brightest examples (Pettini et al. 1998, 2001; 
ies: abundances - galaxies: kinematics - infrared: galaxies Kobulnicky & Koo 2000; Teplitz et al. 2000; Shapley et al. 



. , 2001). Additionally, the 1.5 < z < 2.5 redshift interval remains 

^ ■ unexplored because of the lack of strong spectral features that 

1 . Introduction can be used to identify such sources with visible spectrographs. 



With the recent advent of near-infrared (NIR) spectrographs IMassive clusters acting as Gravitational Telescopes (GTs) 

on 8- 10m class telescopes, a quantitative study of the rest- constitute a powerful tool in the study of high-redshift galax- 

frame optical properties of high-redshift galaxies is now possi- ies. They have been successfully used over a wide range of 

ble. Indeed, NIR spectroscopy allows one to compare directly wavelengths, ranging from the UV to the sub-mm (e.g. Small 

the physical properties of high-redshift galaxies (star forma- et al. 1997; Bezecourt et al. 1999; Ebbels et al. 1998; Pello et 

tion rate, reddening, metallicity, kinematics, virial mass) with al. 1999; Altieri et al. 1999; Ivison et al. 2000) and allowed 

galaxies in the local Universe through line-based indicators. recently to detect the most distant galaxies (z > 5; Ellis et 

The pioneering work by Pettini et al. (1998, 2001) has al. 2001; Hu et al. 2002). The large magnification factors of 

shown that the rest-frame optical properties of the brightest galaxies that are close to the critical lines (typically 1 to 3 mag- 

Lyman break galaxies (LBGs) at z ^ 3 are relatively uniform, nitudes) can be used to probe the physical properties of intrin- 

However, the sample of LBGs observed in the NIR is small sically faint high-redshift galaxies, which would otherwise be 

beyond the limits of conventional spectroscopy (e.g. Ebbels et 

Send offprint requests to: M. Lemoine-Busserolle; ^1. 1996; Pello et al. 1999, Mehlert et al. 2001). 
marie . lemoine@ast . obs-mip . f r 

* Based on observations collected at the Very Large Telescope Large samples of 2;>2.5 star-forming galaxies have be- 

(Antu/UTl), European Southern Observatory, Paranal, Chile (ESO come available in the recent years, mainly through the Lyman 

Programs 64.O-0439, 67.A-0466) break technique (Steidel et al. 1996, 1999). Our sample of 
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highly-magnified, high-redshift galaxies extends these sam- 
ples towards intrinsically fainter galaxies, allowing us to probe 
the physical properties of galaxies that are 1 to 3 magnitudes 
fainter than the present LBGs studies (Steidel et al. 1996, 1999; 
Pettini et al. 1998, 2001). The sample has been selected from 
lensing clusters with well constrained mass distributions, first 
through photometric redshift and lens-inversion techniques, 
and later confirmed spectroscopically (e.g. Kneib et al. 1996; 
Ebbels et al. 1998; Campusano et al. 2001). Lensed galaxies 
for spectroscopic follow-up are chosen to be close to the high- 
redshift critical lines in order to obtain the largest magnification 
(typically ~ 2 magnitudes). 

In this paper, we present NIR spectroscopy on two high- 
redshift galaxies (z ~ 1.9) that are lensed by the massive clus- 
ter AC 114 (z = 0.312). AC 114 is a rich source of multiply- 
imaged, high-redshift galaxies (Small et al. 1995; Natarajan et 
al. 1998, Campusano et al. 2001). The two lensed galaxies, 
named S2 and A2 according to Natarajan et al. 1998, corre- 
spond to multiple images of two different sources at very sim- 
ilar redshifts. A mass-model, first obtained by Natarajan et al. 
(1998) and later improved by Campusano et al. (2001), allows 
us to recover the intrinsic properties of S2 and A2. Thanks to 
the large magnification of these sources, we have measured for 
the first time several important emission lines from NIR spec- 
troscopy of high-redshift galaxies fainter than Mb = — 21. 

The plan of the paper is as follows. In Sections ^ and |[ 
we summarize the observations and describe the steps in con- 
verting the two dimensional raw data into calibrated one di- 
mensional spectra. Emission-line intensities and FWHMs are 
presented in Section ^ and analysed in Section |[ Star forma- 
tion rates, 0/H and N/0 abundance ratios and kinematics are 
derived here. In Section ^ these new data are combined with 
existing samples of LBGs at z ~ 3 in order to understand 
where high-redshift galaxies lie in two fundamental scaling re- 
lations: N/0 versus 0/H abundance ratios and the metallicity- 
luminosity relation. Finally, the main results of this work are 
summarized in Section 0. Throughout this paper, we assume 
a cosmology with Qq = 0.3, A = 0.7 and Hq = 70 km s~^ 
Mpc~^ 

2. Observations 

The data were obtained with the NIR spectrograph ISAAC, 
which is located on the Nasmyth B focus of Antu (VLT-UTl), 
during UT 2001 September 26-27 (period 67). The data were 
obtained with the Short- wavelength channel (Cuby et al. 1999), 
which uses a 1024 x 1024 Hawaii Rockwell array, and the 
medium resolution grating. The slit width was V\ 

Because the redshifts and spectral energy distributions 
(hereafter SED) of the sources in the cluster core were known 
before the run, we could optimize the use of the grating set- 
tings to search for the most relevant spectral features. To cover 
Ha, H/3, [Oiii] AA4959,5007, and [On] AA3726,3728, three 
settings with the ISAAC medium resolution grating were used. 

The ISAAC 120'' long slit allowed us to aUgn at least 2 ob- 
jects per setting in the inner region of the cluster. In the present 
case, we targeted the two gravitationally amplified sources 
AC114-S2 and AC114-A2, both of them at redshift z - 1.9, 




• 



http : / /www . hq . eso . org/ instrument s/isaac/ index . html 



Fig. 1. HSTAVFPC2 image (i?-band/F702W) of the core of 
AC 114. The position of the ISAAC slit used in our observa- 
tions is displayed. Circles mark the position of the high redshift 
galaxies ACl 14-S2 and ACl 14-A2. 

Table 1. Summary of ISAAC Observations 





[Oil] 


[OIII] 


Ha 


Atm. window 


Y 


H 


K 


A (A) observed 


10500-10840 


13900-14500 


18400-19400 


A (A) rest-frame 


3662-3780 


4848-5057 


6417-6766 


Exposure (s) 


8 X 900 


8 X 900 


8 X 900 



together with the offset star. Figure |l] shows an HSTAVFPC2 
image (F702W) of ACl 14 and the long-slit configuration. 

Table [l| summarizes the wavelength regions that were 
selected to contain the redshifted emission lines of [Oil] 
AA3726,3728, [Oiii] AA4959,5007, Ha, and [Nil] A6584. The 
integration times and the airmasses are also listed. 

The ISAAC observations were performed in beam- 
switching mode. The objects were observed at two slit posi- 
tions, which we label A and B. After an A-B-B-A sequence, 
the object was re-acquired at a different slit position, and the 
sequence was repeated. The seeing varied between ~ 0.5 and 
~ 1.5 arcsec and airmass varied between 1.2 and 1.4. 

3. Data reduction and calibration 

We used IRAF with reference to the ISAAC Data Reduction 
Guide]] to convert the two-dimensional raw data into calibrated 
one-dimensional spectra. The data are first split into AB pairs, 
and the first step in the reduction process (the first sky subtrac- 
tion) is to subtract one frame from another in each pair. This 
step removes both the zero level offset of the array and the sky 

^ ISAAC Data Reduction Guide: 
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Fig. 2. Two-dimensional images of the ISAAC spectra of AC114-S2 in the F (1.0 to 1.1 jam), H, and K bandpasses. Positions 
of the nebular emission lines are labeled. 



emission, and results in images with two spectra, one positive 
and one negative. After trimming and flat fielding, we applied a 
two-dimensional wavelength calibration. The wavelength cal- 
ibration, which was derived from the OH lines (Rousselot et 
al. 2000), also corrects for slit curvature. The next step, the 
second sky subtraction, removes any residual sky features re- 
maining after the first sky subtraction. During this process, 
each image is multiplied by —1 and added back to itself af- 
ter a suitable shift. It results in images that have one positive 
spectrum and two negative spectra either side of the positive 
spectrum. The resulting spectra, one for each AB pair, are then 
combined. As an example, the two dimensional spectrum of S2 
is shown in Figure |[ The emission lines [Oil] AA3726,3728, 
[O III] AA4959,5007, and Ha are all clearly detected. 

We extracted a one dimensional spectrum from two- 
dimensional data without fitting the sky since it would only add 
noise. Removing telluric lines is done by dividing the object 
spectrum with that of a telluric standard. The telluric standard 
was observed with the same instrument setup, immediately af- 
ter the science targets, and it was reduced in the same way 
and with the same calibration frames as for the science targets. 
Since the telluric standards were hot stars, a blackbody curve 
was used to model the continuum of these standards. The black- 
body temperature was fixed according to the spectral type of the 
star. The accuracy of telluric line removals has been checked by 
computing the observed [Oiii] A5007/[Oiii] A4959 line ratio 
(see Table g). For S2, we find [Oiii] A5007/[Oiii] A4959 = 
2.5 ± 0.4 which, given the measurement uncertainties, agrees 
with theoretical value of 2.86 (e.g. Osterbrock 1989). Figures ^ 



and y show the extracted one-dimensional infrared Y, H, and 
K band spectra of S2 and A2 respectively. 

4. Emission line measurements 

The emission lines were fitted with one dimensional Gaussians. 
The integrated fluxes (or 2a upper limits when appropriate) 
and velocity widths are listed in Table ^. The intrinsic velocity 
width of emission lines was corrected for the instrumental pro- 
file, which is determined from the width of night-sky lines. The 
redshifts measured for S2 and A2 are 1.867 and 1.869 respec- 
tively. In the case of S2, this value agrees with previous mea- 
surements (Small et al. 1995; Campusano et al. 2001). In the 
case of A2, there is disagreement between the present secure 
value and the estimate by Campusano et al. (2001), z = 1.687. 
Taking into account the correction due to gravitational lens- 
ing, the distance between the two objects is about 4" on the 
source plane at z = 1.869, corresponding to ~ 34 kpc with 
the adopted cosmology. As they have similar redshifts, S2 and 
A2 are thus very close in the source plane and are probably as- 
sociated with one another. Given the redshifts of A2 and S2, 
11(3 is not detected because it is heavily obscured by telluric 
absorption between the Y and H bands. 

The lensing model by Campusano et al. (2001) was revised 
to take into account the correct value for the redshift of A2. 
The resulting magnification factors, 2.0 magnitudes for S2 and 
1.7 magnitudes for A2, are similar to the previous estimates 
(Campusano et al. 2001) and within the model errors (~ 0.1 
magnitudes). 
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Table 2. Spectral Measurements: Line fluxes and FWHM 



Emission Line 




AC114-S2 


AC114-A2 


Comments 


[0 II] A3726 : 










F(xlO~-^^ ergs s~-^ cm" 




L8±0.5 




(1) 


FWHM (A) 




5.4 ±0.7 




(3) 


cr = FWHM/2.35 (km s" 


1) 


65 ±8 




(3) 


[0 II] A3728 : 






[On] AA3726,3728 : 




F(xlO~-^^ ergs s~-^ cm' 




L9±0.4 


L5±0.4 


(1) 


FWHM (A) 




4.2 ±0.4 


11.9 ± 1.8 


(3) 


a = FWHM/2.35 (km s" 




50 ±5 


142 ± 22 


(3) 


[Olll] A4959 : 










F(xlO~^'^ ergs s"^ cm" 




3.4 ±0.4 


< 0.9 (2) 


(1) 


FWHM (A) 




6.9 ±0.5 




(4) 


cr = FWHM/2.35 (km s" 




62 ±5 




(4) 


[Oiii] A5007 : 










F(xlO~-^'^ ergs s~^ cm" 




8.6 ±0.4 


2.5 ±0.4 


(1) 


FWHM (A) 




6.4 ±0.2 


9.5 ±0.4 


(4) 


a = FWHM/2.35 (km s" 




57 ±2 


84 ±3 


(4) 


Ha 










F(xlO~"^^ ergs s~-^ cm" 




15.6 ±0.9 


7.7 ±0.8 


(1) 


FWHM (A) 




6.6 ±0.5 


18.5 ±0.9 


(5) 


cr = FWHM/2.35 (km s" 




45 ±3 


125 ±6 


(5) 


[Nll] A6584 : 










F(xlO~-^^ ergs s~-^ cm" 




< 0.8 (6) 


3.3 ±0.4 


(1) 


FWHM (A) 






13.6 ±0.9 


(5) 


cr = FWHM/2.35 (km s" 






92 ±6 


(5) 



COMMENTS: (1) Measurement by Gaussian fit, corrected for the gravitational magnification factor of 6.31 for ACl 14-S2 and 4.79 for ACl 14- 
A2. (2) Assuming F5007/2.9. (3) Intrinsic FWHM after correction for the instrumental profile, modeled as a Gaussian with FWHM=2.3A. (4) 
Intrinsic FWHM after correction for the instrumental profile, modeled as a Gaussian with FWHM=4.3A. (5) Intrinsic FWHM after correction 
for the instrumental profile, modeled as a Gaussian with FWHM=6.8A. (6) 2a upper limit. 



5. Analysis 

5. 1. Star formation rates 

We have compared the star formation rate (SFR) obtained from 
the UV continuum flux with the value derived from the Ha 
emission line, which is the indicator used for galaxies in the lo- 
cal Universe (e.g. Gallego et al. 1995). The rest-frame 1500A 
flux for S2 and A2 has been derived from the UV rest-frame 
FORSl spectra obtained by Campusano et al. (2001). The spa- 
tial sampling of these spectra is similar to the ISAAC spec- 
tra, thus corresponding to the same physical region in the two 
cases. We have checked these values independently using the 
best-fit SEDs derived from U BVRIJK photometry (see de- 
tails in Campusano et al. 2001). This step was important for 
A2, because the FORSl spectra have a relatively low S/N ra- 
tio. Star formation rates have been computed via the following 
equations (Kennicutt 1998): 

SFRuv(M0yr-i) = 1.4x 10-2%50o(ergs Hz-^), (1) 

SFRHc.(M0yr-i) = 7.9 x Lna{eigs s-^) (2) 

For our adopted cosmology, the SFRs deduced from the Ha 
luminosities are 30.3 and 15.0 M© yr"^ for AC114-S2 and 
AC114-A2, respectively. These values are 7-15 times higher 



than those deduced from the continuum flux at 1500A (4.4 
yr-i for AC114-S2and 1.1 M© yr-^ for AC114-A2) when no 
correction for extinction is applied. If we require that SFRhck 
= SFRuv, then the extinction coefficients are E(B-V)=0.3 for 
AC114-S2 and E(B-V)=0.4 for AC114-A2. With these values, 
the corrected SFRs become 73 M© yr"^ for AC114-S2 and 
51 M0 yr~^ for AC114-A2, applying the reddening law of 
Calzetti et al. (2000). The reddening value of S2 is in good 
agreement with that derived from the UV rest-frame spectrum 
(Le Borgne et al., in preparation). 

These extinctions coefficients are larger than the typical 
value for LBGs up to z ~ 3 (Steidel et al. 1999), and corre- 
spond to the upper envelope of values derived for these objects 
(Shapley et al. 2001). In the remainder of this paper, all line 
fluxes, and hence all the quantities that are derived from these 
fluxes, will be corrected for extinction. 

5.2. Chemical abundances 

We have calculated oxygen-to-hydrogen (0/H) and nitrogen- 
to-oxygen (N/0) abundance ratios for each galaxy using the 
measured emission-line fluxes reported in Table ^. Kobulnicky 
et al. (1999) have shown that emission-line ratios integrated 
over the whole galaxy can provide a reliable indication of the 
oxygen and nitrogen abundances in high-redshift galaxies. 



M. Lemoine-Busserolle et al.: Physical properties of two low-luminosity z ^ 1.9 galaxies 



5 



Table 3. Results of the chemical abundance analysis. 0/H and N/0 abundance ratios using different calibrations for the galaxies 
S2 and A2 in the lensing cluster AC 114. Two cases are considered: without extinction correction and assuming a reddening 
coefficient E(B-V)=0.3 and 0.4 for S2 and A2 respectively. Preferred 0/H abundance ratios are indicated in boldface. 







AC114-S2 




AC114-A2 


51 li 1^1*51 ti mi 

y^aLl UiClLlUll 


E(B-V) 




0.0 


0.3 


U.U 






12 + log(0/H) 














lower branch 




7.17 


7.42 


6.88 


7.25 


from R23 (Pilyugin 2000) 


lower branch 




7.34 


7.66 


7.16 


7.56 


from R23 and O32 (Kobulnicky et al. 1999) 


upper branch 




8.81 


8.63 


8.89 


8.68 


from P and R23 (Pilyugin 2001) 


upper branch 




8.91 


8.78 


8.99 


8.86 


from R23 and O32 (Kobulnicky et al. 1999) 


Secondary abundance indicators 










log([Nll]/Ha) 


< 


-1.29 < - 


-1.29 


-0.36 


-0.36 




log([Nll]/Oll]) 


< 


-0.66 < - 


-0.97 


0.36 


-0.05 




log(N/0) 


< 


-1.32 < - 


-1.51 


-0.42 


-0.68 


from Thurston et al. (1996) 



AC114-S2 ISAAC SZ filter 

[0II]A3726 A3728 



-A2 ISAAC SZ fi 




1.420x10^ 1.430x10^ 

AC1 14-S2 ISAAC SK f 




1.850x10^ 1.860x10' 




Fig. 3. The extracted one-dimensional Y, H, and K band spec- Fig. 4. The extracted one-dimensional Y, H, and K band spec- 
tra of AC114— S2. The dashed line shows the 1 a error spec- tra of AC114— A2. The dashed line shows the 1 a error spec- 
trum. Spectra have been smoothed to the instrumental resolu- trum. Spectra have been smoothed to the instrumental resolu- 



tion. 



tion. 



5.2.1. Oxygen abundance 

Emission lines are the primary source of information regarding 
chemical abundances within Hii regions. The "direct" method 
(referred as the Tq method) for determining chemical com- 
positions requires the electron temperature and the density of 
the emitting gas (e.g., Osterbrock 1989). Unfortunately, a di- 
rect heavy element abundance determination, based on mea- 
surements of the electron temperature and density, cannot be 
obtained for distant galaxies. The [O ill] A4363 auroral line. 



which is the most commonly applied temperature indicator in 
extragalactic Hii regions, is typically very weak and rapidly 
decreases in strength with increasing abundance. 

Given the absence of [O ill] A4363 in our faint spectra, al- 
ternative methods for deriving nebular abundances that rely on 
the bright lines alone must be employed. Empirical methods 
to derive the oxygen abundance exploit the relationship be- 
tween 0/H and the intensity of the strong lines via the param- 
eter R23 = ([Oil] A3726+[Oii] A3728-F[Oiii] A4959+[Oiii] 
A5007)/H/5 both for metal-poor (Pagel, Edmunds, & Smith 
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1980; Skillman 1989; Pilyugin 2000) and metal-rich (Pagel et 
al. 1979; Edmunds & Pagel 1984, McCall et al. 1985; Dopita 
& Evans 1986; Pilyugin 2001) Hll regions. 

The traditional R23 method typically results in metallicities 
within 0.2 dex of those calculated with the more accurate Tg 
abundances. Pilyugin (2000, 2001) improves on the standard 
R23 method by introducing a P factor, which is calculated from 
the strong oxygen lines. The P factor replaces the temperature 
as a descriptor of the conditions in the nebula and results in 
metallicities within 0.1 dex of those calculated with Tq method. 

We thus estimates the oxygen abundance of the two galax- 
ies from the following Pilyugin 's calibrations: 



12 + \og{0/H)R,, = 6.53 + 1.401og(i?23) 



(3) 



12 + \og{0/H)p 



R23 + 54.2 + 59.45P + T.SIP^ 



6.07 + 6.71P + 0.371P2 



0.243i?23 
(4) 



dance in the ionized gas. In this regime, the ionization parame- 
ter also becomes important (e.g., McGaugh 1991). 

Several methods have been proposed to break this degen- 
eracy. Contini et al. (2002) showed that the best abundance 
indicators are the [Nil] A6584/Ha and [Nil] A6584/[Oii] 
AA3726,3728 line ratios which are virtually independent of 
the ionization parameter. Both line ratios can thus be used si- 
multaneously to discriminate between the lower and the upper 
branches in the 0/H vs. R23 relation. The abundance indica- 
tors derived for S2 and A2 are listed in Table ||. They indicate 
clearly that S2 lies on the lower branch of the 0/H vs. R23 
relationship, whereas A2 lies on the upper branch. Additional 
evidence for the low metallicity in S 2 is given by the strength of 
the C IVA1550 absorption line measured in the rest- frame UV 
FORSl spectrum of this object (Le Borgne et al., in prepara- 
tion). Using the relation between 12 + log(0/H) and EW(C IV) 
derived by Mehlert et al. (2002) using Heckman et al. (1998) 
data, we deduce a metallicity Z ~ 0.25 Z©, which is com- 
patible with the value estimated using the strong optical line 
method. 



where P = ([Oiii] A4959+[Oiii] A5007)/([Oii] A3726+[Oii] 
A3728+[Oiii] A4959+[Oiii] A5007) is the excitation parame- 
ter. Equation || refers to the low-metallicity regime and equa- 
tion Q to the high-metallicity one. In order to compare our 
measurements with oxygen abundances reported in the liter- 
ature for nearby and high-redshift star-forming galaxies, we 
also used McGaugh (1991) calibrations with the analytic ex- 
pressions given in Kobulnicky et al. (1999). Oxygen abun- 
dances determined with these calibrations are listed in Table ||. 
Since Ha is detected in S2 and A2 but H/3 is not, we adopt 
= assuming theoretical case B hydrogen re- 

combination ratios for ionized gas with electron temperature 
of 10^ K (Osterbrock 1989). This approximation is very in- 
sensitive to the actual electron temperature but ignores effects 
of dust extinction. To evaluate this effect, we applied the ex- 
tinction coeffici ents E(B-V)=0.3 and 0.4 for S2 and A2 respec- 
tively (see sect. 5J_) and corrected emission-line fluxes. As ex- 
pected, the correction for extinction raises R23 by increasing 
the strength of [O ii] AA3726,3728 relative to H/3. The value 
of 12 + log(0/H) is thus higher (by 0.2 - 0.3 dex) on the 
low-metallicity branch and lower (by 0.1 — 0.3 dex) on the 
high-metallicity branch, depending of the calibration used (see 
Table 1). 

A complication with the "strong lines" method is that the 
dependence of metallicity on R23 is double valued. In the most 
metal-rich Hll regions, R23 is small because metals permit effi- 
cient cooling, reducing the electronic temperature and the level 
of collisional excitation. On the upper, metal-rich branch of the 
relationship, R23 increases as metallicity decreases via reduced 
cooling and elevates the electronic temperature and the degree 
of collisional excitation. However, the relation between R23 
and 0/H changes when 12 + log(0/H) - 8.4 (Z - 0.3 Z©). 
As metallicity decreases below 12 + log(0/H) ~ 8.2, R23 de- 
creases once again. Even though the reduced metal abundance 
further inhibits cooling and raises the electron temperature on 
this lower, metal-poor branch, the intensity of the oxygen emis- 
sion lines drops because of the greatly reduced oxygen abun- 



5.2.2. Nitrogen Abundance 

Nitrogen-to-oxygen abundance ratios (N/0) may be deter- 
mined in the absence of a measurement of the temperature- 
sensitive [Olll] A4363 emission line using the algorithm 
proposed by Thurston et al. (1996). Again, this only re- 
quires the bright [Nil] A6584, [On] AA3726,3728 and [Oiii] 
AA4959,5007 emission lines. The relationship, which is based 
on the same premise as that used in the relationship between the 
oxygen abundance and R23, is calibrated using photoionization 
models. 

First, an estimate of the temperature in the [N ll] emission 
region (t[Nii]) is given by the empirical calibration between 
^[Nii] andi?23 (Thurston et al. 1996): 

t[N II] = 6065+1600(logi?23)+1878(logi?23)^+2803(logi?23)^ 

(5) 

The [Nil] temperature determined from the R23 relation 
can thus be used together with the observed strengths of [N ll] 
A6584 and [On] AA3726,3728 to determine the ionic abun- 
dance ratio N+/0+. Pagel et al. (1992) gave the following for- 
mula based on a five-level atom calculation: 

, N+ , [NII1A6584 ^ ^ 0.726 

(6) 

where the [Nil] temperature is expressed in units of 10^ K. 
Finally, we assume that N/0 = N+/0+ . There has been some 
discussion regarding the accuracy of this assumption (e.g., 
Vila-Costas & Edmunds 1993) but Thurston et al. (1996) found 
that this equivalence only introduces small uncertainties in de- 
riving N/0. The values of N/0 derived for S2 and A2 are listed 
in Table |[ without extinction correction and assuming a red- 
dening E(B-V)=0.3 and 0.4 for S2 and A2 respectively. An up- 
per limit on N/0 is given for S2 because the [Nil] A6584 line 
is not detected. As expected, the extinction correction reduces 
(by ~ 0.2 dex) the N/0 abundance ratio due to the increase of 
[On] AA3726,3728 relative to [Nil] A6584. 
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Table 4. Velocity Dispersions and Virial Masses 



Galaxy a rfj~2 ^1/2 Mvir 

[kms~^] ["] [kpc] [10^° Mq] 

AC114-S2 55 ±6 0.18 1.52 0.53 ± 0.12 

AC114-A2 105 ± 15 0.22 1.85 2.36 ± 0.67 




" From HST WFPC2 images. 

^ = 0.3; A = 0.7; ifo = 70 km s"^ Mpc"^ 

5.3. Kinematics 

We have measured the line-of- sight velocity dispersion a 
of the ionized gas from the brightest emission lines ([Olll] 
AA4959,5007 and Ha) that are corrected for instrumental 
broadening. The values obtained are 55 ± 6 km s~^ for the 
brightest node of ACl 14-S2 and 105 ± 15 km s"^ for ACl 14- 
A2. In order to derive masses for these objects and to com- 
pare with the typical values obtained for LBGs, we estimate 
the virial masses from the following equation: 

Mvir(M0) = 1.16 X lOWi/2 (7) 

where ri/2 is the half-light radius in kpc and a is the velocity 
dispersion in km s~^. This expression corresponds to the ideal 
case of an homogeneous sphere with a uniform density dis- 
tribution. Table ^ summarizes the results. In the present case, 
we have used deep HST-WFPC2 images of this field, obtained 
in the red filter F702W, to compute the half-light radius, after 
correcting for the tangential stretching induced by gravitational 
lensing. In both cases, we have obtained ri/2 ^0.2^'. 

Thanks to the gravitational magnification, which strengths 
the image along the tangential direction, the line profiles of 
ACl 14-S2 are spatially resolved. Figure || displays the raw and 
lens-corrected velocity profiles measured from Ra and [O III] 
A5007 emission lines, after correction for the geometrical po- 
sition of the object across the slit. The detailed lens model has 
been used to derive the tangential magnification along the ob- 
ject. The spatial resolution of the 2D spectra ranges between 
0.5^' and I" without the tangential stretching induced by lens- 
ing (4 to 8 kpc with this cosmology). The velocity gradient 
seems to reach a "plateau" at ~ ±1 kpc from the center of 
the object. In the hypothesis that this velocity gradient traces 
a rotation curve of ±240 (±30) km s~^ in a disk-like ge- 
ometry, the dynamical mass within the inner 1 kpc radius is 
Mdyn = RvyC = 1.3 ± 0.3 X 10^° M0. This value is a factor 
of 2 larger than the mass derived from the velocity dispersion. 

The virial mass derived for ACl 14-A2, and the largest de- 
termination for ACl 14-S2, a few 10^^ Mq, is of the same order 
as the typical values found for LBGs (e.g. Pettini et al. 2001). 
In contrast, AC114-S2 is about one magnitude brighter than 
ACl 14-A2, and one magnitude fainter than the typical LBGs in 
the rest-frame 5-band. Thus, the M/Lb ratio of AC114-S2 in 
solar units ranges between 0.20 and 0.5, the lower value being 
quite similar to that of the LBGs at 2; ~ 3, whereas ACl 14-A2 
has a M/Lb = 2, a ratio ~ 20 times higher than for LBGs. 
After dust extinction correction, we obtain M/Lb ~ 0.15 and 
0.6 for AC114-S2 and AC114-A2 respectively, using the dy- 
namical mass for S2. It is probable that both velocities and 
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Fig. 5. Spatially resolved emission lines in AC114-S2. Top 
panel: HSTAVFPC2 (i?-band/F702W) close image of the 
lensed galaxy ACl 14-S2. Mid panel: Raw velocity profiles (in 
km s~^) versus angular distance (in arcsec), as measured on the 
2D spectra from the central velocities of Ha (open dots) and 
[Olll] A5007 (full dots) emission lines. Bottom panel: Lens 
corrected velocity profiles. 

masses have been underestimated because of the limited sen- 
sitivity of the observations, since we see only the inner cores of 
the galaxies where star formation is the strongest. 

6. Discussion 

Table | summarizes the physical properties for high-redshift 
{z > 2) galaxies. To date there is only one LBG in which 
chemical abundances have been determined with some degree 
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Table 5. Physical properties of high-redshift (z > 2) galaxies. Abundance ratios and SFRs are without extinction correction. 



Galaxy z Mb 12+log(0/H) log(N/0) SFRuy SFRhc Myir (10'^ Mq) 

AC114-S2 1.867 -20.60 7.25 =b 0.2 < -1.32 4.4 30 0.53 ± 0.12 

AC114-A2 1.869 -19.80 8.94 ii= 0.2 -0.42 ± 0.2 1.1 15 2.36 =b 0.67 

CFg 2.313 -22.28 7.8 - 8.7 Z TI 54 T6 

MS 1512-CB58 ... 2.729 -22.04 8.39 ± 0.2 -1.24 ±0.2 16 18 1.8 

Lynx 2-9691 2.888 -23.19 8.3 - 8.8 ... 11 80 5.2 

Q0201-hll3C6 ... 3.055 -22.20 7.6 - 8.8 ... 27 17 1.2 

SSA22aD3 3.069 -22.59 8.0 - 8.6 ... 38 24 

Q1422-h231 D81 . 3.104 -22.91 7.7 - 8.7 ... 45 74 

DSF 2237-^116 C2 3.333 -23.97 7.6 - 8.8 ... 44 75 5.5 

B2 0902-^342 C12 3.387 -22.84 7.6 - 8.8 ... 42 61 



of confidence, the gravitationally lensed galaxy MS 1512-cB58 
at z = 2.729 (Teplitz et al. 2000; Pettini et al. 2002a). Here 
we provide accurate chemical abundance measurements for the 
two high-redshift galaxies: ACl 14-S2 and ACl 14-A2. 

In terms of SFRs without extinction correction, the 
two lensed low-luminosity galaxies A2 and 82 have 
SFRisoo < SFRhc, while in LBGs at z - 3 SFR1500 ^ 
SFRhq (within 15%). Note that a factor of two still remains 
between the two SFRs of S2 when using the latest calibrations 
by Rosa-Gonzalez et al. (2002). As shown in Section 5J_, a 



large E(B-V) value is needed to reconcile the UV/optical SFR 
values for these two objects. Unfortunately, we cannot derive 
an independent reddening estimate for S2 and A2, e.g. using 
the Balmer decrement Ha/H/3.. 

6. 1. The metallicity-luminosity relationship 

We can study the evolution with redshift of the fundamen- 
tal scaling relation between galaxy luminosity and metallicity 
(Contini et al. 2002 and references therein). For nearby galax- 
ies, this relation extends over ~ 10 magnitudes in luminosity 
and ~ 2 dex in metallicity. 

The high-redshift sample is best compared to nearby galax- 
ies where metallicities are derived using the same empir- 
ical strong line method. Figure ^ shows the relations for 
nearby i) "normal" irregular and spiral galaxies (solid line; 
Kobulnicky & Zaritsky 1999) and ii) starburst galaxies (dashed 
fine; Mouhcine & Contini 2002; Melbourne & Salzer 2002), as 
well as the high-redshift sample from Table ^ Three samples of 
intermediate-redshift galaxies are also shown for comparison: 
UV-selected star-forming galaxies at z ~ 0.1 — 0.4, emission- 
line galaxies at z~ 0.1 — 0.5 (Kobulnicky & Zaritsky 1999), 
and luminous compact emission-line galaxies atz^O.S — 0.7 
(Hammer et al. 2001). 

It is immediately obvious that high-redshift galaxies do 
not conform to today's luminosity-metallicity relation for both 
"normal" and starburst galaxies. Even allowing for the un- 
certainties in the determination of 0/H, high-redshift galaxies 
have much lower oxygen abundances than one would expect 
from their luminosities. This result, already revealed by pre- 
vious studies (Kobulnicky & Koo 2000; Pettini et al. 2001; 
Contini et al. 2002), is secured with the addition of the low- 
luminosity and low-metallicity galaxy S2. 



One interpretation of this result is that high-redshift galax- 
ies are undergoing strong bursts of star formation which raise 
their luminosities above those of nearby galaxies with simi- 
lar chemical composition. Another possibility is that the whole 
metallicity-luminosity relation is displaced to lower abun- 
dances at high redshifts, when the Universe was younger and 
the time available for the accumulation of the products of 
stellar nucleosynthesis was shorter. It should be possible to 
quantify this effect by measuring metallicities in samples of 
galaxies at different redshifts (e.g. Mehlert et al. 2002). As 
shown in Figure |^, intermediate-redshift galaxies seem to fol- 
low the luminosity-metallicity relation derived for nearby star- 
burst galaxies (dashed line). Unfortunately, to date there is no 
sample of galaxies with known chemical abundances to fill the 
gap between z ~ 0.5 and z ^ 2. 

The location of A2 is more surprising. It does not follow the 
trend of high-redshift objects. Instead, it lies on the luminosity- 
metallicity relation of nearby objects and has the highest metal- 
licity of any high-redshift (z > 2) galaxy. 

We compute the intrinsic Mi from the observed i^T-band 
magnitudes (2.2 fim corresponds to 7600 A at z ~ 1.9) for S2 
and A2. We found the same value Mi = —23.2 for S2 and A2. 
This suggests that S2 and A2 may have the same stellar masses, 
since in the / band we have access to the old stellar population 
which dominates the stellar mass. But the relatively high nitro- 
gen abundance in AC114-A2 and MS 1512-cB58 suggests the 
presence of an older population of stars than in ACl 14-S2. 

In Figure we plot the oxygen abundance (12 + log(0/H)) 
against the line-of- sight velocity dispersion (a) for S2 and A2 
and other relevant samples of nearby (Telles & Terlevich 1997; 
Richer & McCall 1995; Zaritsky, Kennicutt & Huchra 1994; 
Oey & Kennicutt 1993), intermediate (Kobulnicky & Zaritsky 
1999; Hammer et al. 2001) and high redshift (see Table |) 
galaxies. This diagram corresponds to a crude mass-metallicity 
sequence. The location of AC114-A2 is fully compatible with 
that of low-redshift galaxies. However, the oxygen abundance 
of ACl 14-S2 is much smaller than the corresponding value for 
low-z galaxies of similar velocity dispersions. The situation is 
even more dramatic if we consider the velocity gradient derived 
under the disk-like rotation curve hypothesis (see sect. 5.3). 
Taken at face value, it would appear that the slope of the mass- 
metallicity relation for galaxies with z > 2 is different from 
that of low- 2; galaxies. 
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Fig. 6. Oxygen abundance (12 + log(0/H)) as a function of 
5-band absolute magnitude. The metallicity-luminosity rela- 
tions for nearby "normal" (solid line; Kobulnicky & Zaritsky 
1999) and starburst (dashed line; Mouhcine & Contini 2002; 
Melbourne & Salzer 2002) galaxies are shown. The location 
of the two z ^ 1.9 galaxies S2 and A2 in the lensing cluster 
AC 1 14 is shown without extinction correction (empty squares) 
and assuming a reddening E(B-V)=0.3 and 0.4 for S2 and A2 
respectively (filled squares). We show for comparison sam- 
ples of intermediate-redshift galaxies: UV- selected galaxies at 
z ~ 0.1 — 0.4 (circles; Contini et al. 2002), emission-line galax- 
ies at z ~ 0.1 — 0.5 (stars; Kobulnicky & Zaritsky 1999), and 
at z ~ 0.5 — 0.7 (crosses; Hammer et al. 2001). The loca- 
tion of high-redshift (z ~ 3) Lyman break galaxies (LBGs) 
is shown as a blue box encompassing the range of 0/H and 
Mb derived for these objects (Pettini et al. 2001). High-redshift 
(1.4 < z < 3.4) galaxies with oxygen abundances and 5-band 
magnitudes derived from rest- frame UV spectra (Mehlert et al. 
2002) are shown as blue squares. 

6.2. The N/O versus O/H relationship 

In Figure ^, we examine the location of S2 and A2 in the N/O 
versus O/H relationship. The behavior of N/O with increasing 
metallicity offers clues about the history of chemical evolution 
of galaxies and the stellar populations responsible for produc- 
ing oxygen and nitrogen. The location of the two z ^ 1.9 
galaxies AC114-S2 and AC114-A2 are shown without extinc- 
tion correction (empty squares) and assuming a reddening E(B- 
V)=0.3 and 0.4 for S2 and A2 respectively (filled squares). For 
comparison, we show in Figure |^ the location of nearby star- 
forming galaxies: starburst nucleus galaxies (SBNGs) and Hii 
galaxies (see Contini et al. 2002 for references). We also plot 
samples of intermediate-redshift galaxies: UV- selected galax- 
ies 3.1 z ~ 0.1 — 0.4 (Contini et al. 2002), and emission-line 
galaxies at 2; ~ 0.1 — 0.5 (Kobulnicky & Zaritsky 1999). The 



Fig. 7. Oxygen abundance (12 + log (O/H)) as a function of the 
line-of-sight velocity dispersion (a). The location of ACl 14-S2 
and AC114-A2 is shown without extinction correction (empty 
squares) and assuming a reddening E(B-V)=0.3 and 0.4 for 
S2 and A2 respectively (filled squares). We show for compari- 
son samples of nearby galaxies: Hii galaxies (triangles; Telles 
& Terlevich 1997), dwarf galaxies (crosses; Richer & McCall 
1995), late-type spirals (squares; Zaritsky, Kennicutt & Huchra 
1994) and early-type spirals (circles; Oey & Kennicutt 1993). 
The dashed line is a linear fit to these samples. Intermediate- 
redshift galaxy samples are: emission-line galaxies at z ~ 
0.1 - 0.5 (stars; Kobulnicky & Zaritsky 1999), and Sit z ^ 
0.5 — 0.7 (crosses; Hammer et al. 2001). The location of high- 
redshift (z ~ 3) Lyman break galaxies (LBGs) is shown as 
a blue box encompassing the range of O/H and a derived for 
these objects (Pettini et al. 2001). 

location of Damped Lya systems (Pettini et al. 2002b) is also 
shown. 

The sample of high-redshift (z > 2) galaxies with mea- 
sured N/O abundances is still very small: A2 and S2 from 
this paper, and MS 1512-cB58 at z - 2.7 (see Table |). 
Surprisingly, these three galaxies have very different locations 
in the N/O vs. O/H diagram. S2 is a low-metallicity object 
(12 + log(0/H) - 7.25; Z - 0.03 Z©) with a low N/O ra- 
tio (N/O < —1.32), similar to those derived in the most metal- 
poor nearby Hll galaxies. In contrast, A2 is a metal-rich galaxy 
(12 + log(0/H) - 8.94; Z - 1.3 Zq) with a high N/O abun- 
dance ratio (N/O ~ —0.42), similar to those derived in the most 
metal-rich massive SBNGs. The position of MS 1512-cB58 is 
intermediate between these two extremes showing abundance 
ratios typical of low-mass SBNGs and intermediate-redshift 
galaxies (see Fig. ^. 

A natural explanation for the variation of N/O at constant 
metallicity is the time delay between the release of oxygen and 
that of nitrogen into the ISM (e.g. Contini et al. 2002, and refer- 
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Fig. 8. Nitrogen-to-oxygen (N/0) abundance ratio as a func- 
tion of oxygen abundance (12 + log(0/H)). The location of 
the two z ^ 1.9 galaxies S2 and A2 in the lensing cluster AC 
1 14 is shown without extinction correction (empty squares) and 
assuming a reddening E(B-V)=0.3 and 0.4 for S2 and A2 re- 
spectively (filled squares). Two comparison samples of nearby 
star-forming galaxies are shown (see Contini et al. 2002 for ref- 
erences): starburst nucleus galaxies (dot - short dash line) and 
Hll galaxies (short dash - long dash line). We show, for com- 
parison, samples of intermediate-redshift galaxies: UV-selected 
galaxies at z ~ 0.1 — 0.4 (circles; Contini et al. 2002), and 
emission-line galaxies at z~ 0.1 — 0.5 (stars; Kobulnicky & 
Zaritsky 1999). The location of Damped Lya systems (trian- 
gles, filled symbols are upper limits; Pettini et al. 2002b) is 
shown for comparison. Theoretical curves for a primary (dot- 
ted line), a secondary (dashed line), and a primary + secondary 
(solid line) production of nitrogen (Vila-Costas & Edmunds 
1993) are shown. 



ences therein), while maintaining a universal IMF and standard 
stellar nucleosynthesis. The "delayed-release" model assumes 
that star formation is an intermittent process in galaxies (e.g. 
Edmunds & Pagel 1978; Garnett 1990; Coziol et al. 1999) and 
predicts that the dispersion in N/0 is due to the delayed release 
of nitrogen produced in low-mass longer-lived stars, compared 
to oxygen produced in massive, short-lived stars. 

Following this hypothesis and new chemical evolution 
models (Mouhcine & Contini 2002), we might interpret the lo- 
cation of S2 and A2 in the N/0 versus 0/H diagram in terms of 
star formation history and evolutionary stage of these galaxies. 
The low 0/H and N/0 abundance ratios found in S2 might in- 
dicate a relatively young age for this object which experienced 
two or three bursts of star formation at most in the recent past. 
A2 seems on the contrary much more evolved. The location 
of this galaxy in the N/0 versus 0/H diagram indicates a rela- 



tively long star formation history with numerous powerful and 
extended starbursts. 

7. Conclusions 

We have presented in this paper the firsts results of our NIR 
spectroscopic survey of highly magnified high-redshift galax- 
ies in the core of lensing clusters. Our first targets were ACl 14- 
S2 and ACl 14-A2, two lensed sources at z ~ 1.9. We have ob- 
tained, for the first time, the rest-frame optical emission lines 
of galaxies 1-2 magnitudes fainter than the LBGs studied in 
field surveys (e.g. Pettini et al. 1998, 2001; Kobulnicky & Koo 
2000). The main results obtained in this paper are the follow- 
ing: 

1 . The SFRs derived from the Ha emission-line are system- 
atically higher than the values deduced from the 1500A 
flux without dust extinction correction. The reddening val- 
ues needed to obtain SFRna/SFRuv ~ 1 using the Calzetti 
et al. (2000) reddening law are E(B-V)=0.29 for ACl 14-S2 
and E(B-V)=0.40 for ACl 14-A2, leading to corrected SFRs 
of 73 M0 yr~^ and 51 Mq yr~^ respectively. These results 
suggest that large dust extinction corrections are needed for 
intrinsically faint galaxies compared to field (brighter) LEG 
samples. 

2. The large spectral coverage obtained, from [On] 
AA3726,3728 to Ha+[Nii] A6584, allows us to set 
strong constraints on the chemical abundance ratios (0/H 
and N/0) of these two faint galaxies. The behavior of S2 
and A2 in terms of metallicity is different, and the two 
objects are also different from typical LBGs at z ~ 3, 
suggesting a different star formation history for galaxies of 
different luminosities. 

3. The virial masses derived from the line-of-sight velocity 
dispersions are of the order of 5 x 10^ Mq for ACl 14- 
S2 and - 10^° M© for AC114-A2. The line profiles of 
AC114-S2 are spatially resolved, leading to a velocity gra- 
dient of ±240 km s~^, which yields a dynamical mass of 
1.3±0.3 X 10^^ M0 within the inner 1 kpc radius, under 
the hypothesis of a rotation curve with disk-like geometry. 
Although the M/Lb ratio of AC114-S2 is quite similar, 
within a factor of 2, to the values found for field (~ 1 mag 
brighter) LBGs, the corresponding value for AC114-A2 is 
about one order of magnitude higher. After dust extinction 
correction, we obtain M/Lb ^ 0.15 and 0.6 for S2 and A2 
respectively. 

In terms of SFRs, mass-to-light and chemical abundance 
ratios, the behaviour of S2 and A2 is rather close to that de- 
scribed by Guzman et al. (1997) for their sample of compact 
emission-line galaxies in the HDF at 0.7 ^ z ^ 1.4. S2 would 
correspond to a young star-forming Hll galaxy, whereas A2 is 
more likely an evolved starburst. 

The results obtained on the physical properties of ACl 14- 
S2 and AC114-A2 suggest that high-2: objects of different lu- 
minosities could have quite different star formation histories. 
However, the number of well observed high redshift objects is 
currently very small and larger samples in both redshift (1.5 
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< z < 6) and luminosity are required. This could be achieved 
with the new generation of multi-object NIR spectrographs for 
the 10m class telescopes, such as KMOS on the VLT or EMIR 
on the GTC. 
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